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The weak localization correction to the conductivity in coupled double layer structures is studied 
both experimentally and theoretically. Statistics of closed paths has been obtained from the analysis 
of magnetic field and temperature dependencies of negative magnetoresistance for magnetic field 
perpendicular and parallel to the structure plane. The comparison of experimental data with results 
of computer simulation of carrier motion over two 2D layers with scattering shows that inter-layers 
transitions play decisive role in the weak localization. 



PACS numbers: 73.20Fz, 73.6fEy 



I. INTRODUCTION 

Transitions between 2D layers is one of fundamental 
features of double layer structures. It changes the quan- 
tum corrections to the conductivity, especially in a mag- 
netic field paralleLto the layers. 

It is well knownQ that the interference of electron waves 
scattered along closed trajectories in opposite directions 
(time-reversed paths) produces a negative quantum cor- 
rection to the conductivity. An external magnetic field 
(B) gives the phase difference between pairs of time- 
reversed paths (p = 27t(BS)/$Oj where $o is the quan- 
tum of magnetic flux, S is the area enclosed, and thus 
destroys the interference and results in negative magne- 
toresistance. 

In case of a single 2D layer the influence of a magnetic 
field is strongly anisotropic because all the paths lie in 
one plane. The magnetoresistance is maximal for B || n, 
where n is the normal to 2D layer. When a magnetic 
field lies in the 2D layer plane, Bin, the scalar prod- 
uct (BS) is zero, i.e. the magnetic field does not destroy 
the interference, and the negative magnetoresistance is 
absent in this magnetic field orientation.H 

In coupled double layer structures, the tunneling be- 
tween layers gives rise to the closed paths where an elec- 
tron moves initially over one layer then over another one 
and returns to the first layer. For this paths the product 
(BS) is non-zero for any magnetic field orientation and 
hence the negative magnetoresistance has to appear for 
B _L n as well. 

The magnetic field dependence of the negative mag- 
netoresistance is determined by the statistics of closed 
paths, namely, by the area distribution function, W(S), 
and area dependence of the average length of closed 
paths, L(S)B@ Just these statistic dependencies have 
been studied in single 2D layer structures by analysis 
of negative magnetoresistance measured at B |j n.aC 

The role of inter-layers transitions in weak localiza- 
tion and negative magnetoresistance for B || n for mul- 



tilayer structures (superlattices) was discussed in Ref. [t]. 
The closely related problem concerning the role of inter- 
subbands transitions in quasi-two dimensional structures 
with several subbands occupied was theoretically studied 
in Ref. H 

In this paper we present the results of investiga- 
tions of the negative magnetoresistance in double layer 
GaAs/InGaAs structure for different magnetic field ori- 
entations. We obtain the area distribution functions and 
area dependencies of the average lengths of the closed 
paths using the approach developed in Refs. These 
functions are compared with those obtained from the 
computer simulation of carrier motion when inter-layers 
transitions are accounted for. Close agreement shows 
that just the inter-layers transitions determine the neg- 
ative magnetoresistance in coupled double layer struc- 
tures. 



II. EXPERIMENTAL RESULTS 

The double well heterostructure GaAs/In^Gai_ x As 
was grown by Metal-Organic Vapor Phase Epitaxy 
(MOVPE) on semi-insulator GaAs substrate. The het- 
erostructure consists of a 0.5 /im-thick undoped GaAs 
epilayer, a Si 5-layer, a 75 A spacer of undoped GaAs, 
a 100 A Ino.08Gao.92As well, a 100 A barrier of undoped 
GaAs, a 100 A Ino.osGao.92As well, a 75 A spacer of 
undoped GaAs, a Si <5-layer and 1000 A cap layer of un- 
doped GaAs. The samples were mesa etched (cut??) 
into standard Hall bridges. The measurements were per- 
formed in the temperature range 1.5 — 4.2 K at low mag- 
netic field up to 0.4 T with discrete 10~ 4 T for two ori- 
entations: the magnetic field was perpendicular (B || z) 
and parallel (B || x) to the structure plane (see insert 
in Fig. |l|). Additional high field measurements were also 
made to characterize the structure. It has been found 
that in the structure investigated the conductivity is de- 
termined by the electrons in the wells. Their densities 
have been determined from the Fourier analysis of the 
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Shubnikov-de Haas oscillations and consist of 4.5 x 10 11 
cm~ 2 and 5.5 x 10 11 cm -2 in different wells. The Hall 
mobility was about fx ~ 4200 cm 2 / (V x sec). 
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FIG. 1. Magnetic field dependencies of Aa /Go for different 
magnetic field orientations for T — 4.2 (a), 1.5 K (b). Sym- 
bols are the experimental data, red curves are the simulation 
results. Blue curves are the results of calculations carried out 
according to Ref. |ll| Insert in (a) shows a system of coordi- 
nates. 

The magnetic field dependencies of in-plane magneto- 
conductance 



Aa(B) = a(B) - a(0) = l/p(B) - l/p(0) 



(1) 



at magnetic field perpendicular (Aa(B z )) and parallel 
(Aa(B x )) to the structure plane are presented in Fig. [|. 
One can see that the negative magnetoresistance is ob- 
served for both magnetic field orientations and, in con- 
trast to the case of single layer structures, the effects 
are comparable in magnitude. Analysis of behaviour 
of the conductivity in a wide range of temperatures 
(1.5 < T < 20 K) and magnetic fields (B < 6 T) shows 
that at B < 0.4 — 0.5 T the main contribution to the 
negative magnetoresistance comes from the interference 
correction. In this case the magnetic field dependence 
of the negative magnetoresisjtaiice is determined by the 
statistics of the closed pathsBEl 

Let us apply the method proposed in Refs. to anal- 
ysis of negative magnetoresistance in the double layer 
structure. Using formalism presented in Section II of Ref. 

one can write the expression for conductivity of dou- 
ble layer structure with identical layers for two magnetic 
field orientations as follows 



<r(Bt) 



= cr + Sa(Bi) 
= cr - 47rZ 2 G 



exp 



dS l IwiS, 
(2-KBiSi 



cos 



(2) 



Here, Go = e 2 /(27r 2 7i), <7o is the classical Drude conduc- 
tivity, I = vpT, l v — vpT v , vp is the Fermi velocity, r 
and t v are the momentum relaxation and phase breaking 
time, respectively. The value of L is the function not only 
of S but l v as well. It was defined in Ref. | by Eq.(6). It 
should be mentioned that Eq. (||) is valid at low enough 
probability of inter-layers transitions. 

Thus, for the magnetic field perpendicular to the struc- 
ture plane (B = (0,0, B z )) the magnetoresistance is 
determined by z-component of S only and for parallel 
magnetic field (B = (B x ,0,0)) it is determined by x- 
component of S. 

One can see from Eq. (§) that the Fourier transform 
of Sa(B)/G 



1 

*o J 



,„ 6a(Bi) {2ttBA 
aBi — - — cos 



^l 2 W(Si) exp 



G V *o 
LiSiY 



(3) 



carriers an information on W{Si) and L(Si). Because 
the value of l v tends to infinity when temperature tends 
to zerojj the extrapolation of <&{Si, l v )-vs-T curve to 
T = gives the value of 4ttZ 2 PF(S 1 ). The ratio 1(5*)/^ 
for given Si can be then obtained as \n(Airl 2 W(Si)) — 
ln($($,L,)). 
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FIG. 2. Area (a) and temperature (b) dependencies of 
the Fourier transforms of Sa(B x ). Curves in (b) show the 
extrapolation of T) to T = 0. 

The value Aa(B) = a(B) - a(0), not WB), is experi- 
mentally measured. It is clear from Eqs. (0) and (g) that 
5(i(B) = <t(0) — (To~ Aa(B). To obtain 8a(B), we assume 
that the Drude conductivity <To is equal to the conductiv- 
ity at T=20 K, when the quantum corrections are small. 
Notice that the final results are not sensitive to the value 
of (To practically. Obtaining of the distribution function 
W(S X ) from the experimental Scr(B z ) dependencies for 
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the structure investigated is illustrated by Fig. ||. In left 
panel the Fourier transforms of Sa(B x ) measured at dif- 
ferent temperatures are presented. Right panel shows 
how the <£>-vs-T data have been extrapolated to T = 0. 
The area distribution function W(S Z ) has been obtained 
from the analysis of Sa(B z ) curves in a similar way. 

The results of data processing described above are pre- 
sented in Fig. |a. As is seen the 4tt1 2 W(S z ) dependence 
is close to (25 2 )" 1 for S ~ (0.3 - 5) x 1CT 10 cm 2 . The 
analogous behaviour of the area distribution function was 
obtained for single 2D layer in Ref. ^. The behaviour of 
W(S X ) significantly differs from that of W(S Z ). In par- 
ticular, the W(S X ) curve shows a much steeper decline 
for S > 0.8 x 10~ 10 cm 2 . Other feature of the statis- 
tics of the closed paths in double layer structure is the 
fact that for given S the values of L(S X ) are significantly 
larger than L(S Z ) (Fig. |b). 
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FIG. 3. The area distribution functions of closed paths (a) 
and the area dependence of the L(S X ) / L(S Z ) ratio at T = 1.5 
K (b). The symbols are the experimental data, curves are 
the results of simulation with t = 0.1, red line in (a) is 
(2S) - ^dependence. 

Qualitatively these peculiarities of the statistics of 
closed paths in double layer structures can be under- 
stood if one considers how trajectories with large enough 
length, L ^> l/t, look. They are isotropically smeared 
over xy-plane for the distance ~ their extended 

area in this plane is s z ~ LI. In j/z-plane they have size 
~ \fL~l in y-direction and Zq (where Zq is inter-layer dis- 
tance) in z-direction. So, the extended area in yz-plane is 
s x ~ Zq^/LI. Thus, closed trajectories have significantly 
larger s z than s x , and the s z /s x ratio increases with in- 
creasing s. It is clear that the behaviour of enclosed areas 
S z , S x is analogous. Therefore, for S x = S z the inequal- 
ity W(S Z ) > W(S X ) is valid. The average length of the 
trajectories L(S X ) therewith is greater than L(S Z ). 

As was shown in Ref. ^ the distribution function of 
closed paths, the area dependence of average length of 



closed paths and weak localization magnetoresistance can 
be obtained by computer simulation of a carrier motion 
over 2D plane. 



III. COMPUTER SIMULATION 

The model double layer system is conceived as two 
identical plains with randomly distributed scattering cen- 
ters with a given total cross-section. Every plane is repre- 
sented as a lattice M x M with lattice parameter a. The 
scatterers are placed in a part of the lattice sites with 
the use of a random number generator. We assume that 
a particle moves with a constant velocity along straight 
lines which happen to be terminated by collisions with 
the scatterers. After every collision the particle has two 
possibilities: it passes from one plane to another with a 
probability t and moves over the second plane or it re- 
mains in the plane with probability (1 —t), changing the 
motion direction only. If the trajectory of the particle 
passes near the start point at the distance less than d/2 
(where d is a prescribed value, which is small enough), it 
is perceived as being closed. The projections of enclosed 
algebraic area is calculated according to 
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3 = 1 



Vj +i + Vj 
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Vn + 2/1 



(Xj+l ~Xj) 
(x N - Xi), 



(4) 



N-l 



Sx = j: 



Vj+i + Vj 



(z j+ i ~Zj) + 



3 = 1 



Vn + Vi 



(z N - Zl), 



(5) 



where N is number of collisions for given trajectory, Xj, 
Dj, Zj stand for coordinates of j-th collision, Zj takes the 
value or Zq. Otherwise the simulation details are anal- 
ogous to them described in Ref. || for the case of single 
2D layer system. 

All the results presented here have been obtained us- 
ing the parameters: lattice dimension is 6800 x 6800; the 
number of starts, I s , is 10 6 ; the total number of scatter- 
ers is about 1.6 x 10 5 ; the scattering cross section is 7; 
d = 1; Zq = 18. The value of mean free path computed 
for such a system is about 43 x a. If we suppose the value 
of a equal to 11 A, this model double layer system corre- 
sponds to the heterostructure investigated. Namely, the 
mean free path is equal to the value of I ~ 480 A, and the 
value of Zq is close to the distance between the centers 
of the quantum wells, 200 A. 

The area distribution functions obtained as the result 
of simulation with different inter-layers transition proba- 
bilities are presented in Fig. [|. Let us discuss at first the 
behaviour of W{S Z ) (Fig. @a). For t = 0, the Att1 2 W(S z ) 
curve corresponds to the area distribution function for 
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single layer. For large S this curve goes close to the in- 
dependence, which corresponds to the ideal 2D system in 
the diffusion regime.Ea The deviation, which is evident for 
S z < 10 3 a 2 , is just due to the transition to the ballistic 
regime. It is obviously that for sufficiently large values 
of t, the probability of return to the start point has to be 
twice smaller than that for t = 0. As is seen from Fig. ||a 
even the value t = 0.1 is large enough in this sense: the 
corresponding curve is close to the (2S) _1 -dependence 
practically in whole area range. This is because the long 
trajectories with large number of passes between lay- 
ers give significant contribution to W(S Z ) starting from 
small areas, S z > 0.1 1 2 . For the intermediate value of t 
(t = 0.002, 0.01) the area distribution function is close to 
the S -1 function for small areas and tends to the (2S') _1 
dependence for large ones. 
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FIG. 4. Area distribution functions W(S Z ) (a) and W(S X ) 
(b) as they have been obtained from the simulation procedure 
with different t values. 

The behaviour of W(S X ) contrasts with that of W(S Z ) 
(Fig. [|b). At small S x , W(S X ) depends only weekly on 
S x , whereas at large S x it decreases sharply when S x in- 
creases. Sensitivity of W(S X ) to inter-layers transition 
probability depends on S x value. For small S x values, 
when the area distribution function is mainly determined 
by short closed paths with small number of inter-layers 
transitions, the value of W(S X ) considerably increases 
with increasing t. For large S x , i.e. for paths with large 
number of intcr-layer transitions, W(S X ) weakly depends 
on the transition probability. 

Let us demonstrate how the magnetoresistance of our 
model 2D system changes with changing of the inter- 
layers transition probability. The theoretical Sa(B) de- 
pendencies have been calculated by summing over the 
contributions of every closed path to-the conductivity in 
accordance with following expressiontl 



where k is the length of i-th closed path. The results 
of calculation are presented in Fig. @, where A<j(S,) = 
Sa(Bi) - Sa(0) is plotted against B/B t , B t = hc/(2el 2 ). 
As is seen the changes in magnetic field dependencies of 
negative magnetoresistance with change of inter-layers 
transition probability reflect the variation of area distri- 
bution functions. Indeed, Aa(B z ) depends on t slightly: 
maximal change is less than two times for decreasing t up 
to zero, whereas Aa(B x ) changes drastically. It decreases 
about hundred times, when the value of t decreases from 
0.1 to 0.002. 



IV. DISCUSSION 



Let us compare the calculated area distributions with 
experimental data. One can see from Figs. ||a and^, that 
the behaviour of calculated and experimental W(S Z ) and 
W(S X ) dependencies is close qualitatively. As mentioned 
above, W{S X ) depends on inter-layers transition proba- 
bility significantly stronger than W(S Z ). Therefore we 
have estimated the transition probability comparing the 
calculated and experimental W(S X ) curves. The most 
accordance has been obtained with t ~ 0.1 (see Fig. ||a). 
As seen from the figure, with this value of t the calculated 
W(S Z ) dependence describes the experimental data well. 
Some quantitative inconsistency, especially for W(S X ), is 
evident in Fig. ||a. We believe, this is result of crudity of 
model used. In particular, we supposed the identity of 
both 2D layers. 
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FIG. 5. Calculated magnetic field dependencies of Act for 
different inter-layers transition probability, I /lip = 0.01. 

Let us turn now to magnetic field dependencies of nega- 
tive magnetoresistance. To calculate Aa(B) 7 in addition 
to the inter-layers transition probability it is necessary 
to know the phase breaking length. Using the value of 
t = 0.1 estimated above, we have found that the best 
agreement between theoretical and experimental Aa(Bt) 
dependencies is obtained with l v ~ 3.4 and 1.4 /mi for 
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T=1.5 and 4.2 K, respectively. The Aa(B z ) and Aa(B x ) 
dependencies calculated with these l v values practically 
coincide with those measured experimentally (see Fig. 
|l|) . It should be noted that these values of l v some differ 
from those obtained by-fitting of the Aa(B z ) curves to 
the Hikami expressionild the fit gives l v ~ 4.8 and 1.7 /im 
for T=1.5 and 4.2 K, respectively. The reason of this dif- 
ference is that the Hikami formula was obtained for single 
2D layer, and it is not suitable for analysis of negative 
magnetoresistance in coupled double layers structures. 

Finally, knowing the values of t and l v we are able to 
compare the calculated and experimental area dependen- 
cies of L(S X ) to L(S Z ) ratio (Fig. ^|b). It is seen that the 
experimental ratio is significantly larger than unity as 
well as calculated one. However, the experimental points 
lie somewhat below than calculated curve. We suppose 
that the main reason of such disagreement is dissimilarity 
of the layers in structure investigated. 

After this paper has been prepared for publication, the 
paper by Raichev and Vasilopoulos on the theory of weak 
localization in double quantum .jaells is appeared in Con- 
densed Matter e-Print archive.til Let us apply this the- 
ory to our case. Using the formulae derived in Ref. [ll], 
we have calculated Aa(Bi) dependencies for our struc- 
ture. These dependencies are represented in Fig. [l] by 
dashed curves. As is clearly seen, theory developed in 
Rcf. ^ describes our experimental results only in low 
magnetic fields. The reason is that the calculations in 
Ref. [ll] were carried out in the framework of diffusion 
approximation. It means that two conditions are met. 
The first condition is t < t ? . For the structure inves- 
tigated t/t v — 0.014 — 0.035 for different temperatures, 
and this condition may be considered as fulfilled. Accord- 
ing to the second condition, the magnetic field has to be 
low enough: B <C B t when B || z, or B <C B t l/Zo when 
B || x. In our case B t ~ 0.14 T, l/Zo ~ 2.5 and hence the 
diffusion approximation is applicable, when B <C 0.14 or 
0.35 T depending on the magnetic field orientation. It is 
in this range of magnetic field that the results of Ref. [ll] 
are close to our experimental data. 

Our calculations are valid beyond the diffusion approx- 
imation and therefore they better describe the experi- 
mental results in whole magnetic field range, where weak 
localization correction to the conductivity is dominant. 



V. CONCLUSION 

We have investigated the negative magnetoresistance 
in double layer heterostructures for different magnetic 



field orientations. The information about statistics of 
closed paths has been extracted from the analysis of tem- 
perature and magnetic field dependencies of conductiv- 
ity. Significant difference in area distribution functions, 
W(S X ), W(S Z ), and in average lengths of closed paths, 
L(S X ), L(S Z ), has been found. In order to interpret the 
experimental results, we have investigated the statistics 
of closed paths and negative magnetoresistance using the 
computer simulation of the carrier motion with scatter- 
ing over two 2D layers. Analysis of experimental and 
theoretical results unambiguously shows that in parallel 
magnetic field the negative magnetoresistance in double 
layer structures is determined by inter-layers transitions. 
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